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The classical Yang-Mills equations are solved perturbatively in covariant gauge for a 
collision of two ultrarelativistic nuclei. The nuclei are taken as ensembles of classical color 
charges on eikonal trajectories. The classical gluon field is computed in coordinate space up 
' to cubic order in the coupling constant g. We construct the Feynman diagrams correspond- 

I ing to this field and show the equivalence of the classical and diagrammatic approaches. An 

■ argument is given which demonstrates that at higher orders in g the classical description 

^\ ' of the process breaks down. As an application, we calculate the energy, number, and mul- 

tiplicity distributions of produced soft gluons and reproduce earlier results by Gunion and 
Bertsch and by Kovner, McLerran, and Weigert. 
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I. INTRODUCTION 



^ ■ Ultrarelativistic heavy-ion collisions {^/s ~ 200 AGeV) at the Relativistic Heavy Ion Collider (RHIC) at 

I Brookhaven aim towards an understanding of the properties of nuclear matter under extreme conditions ^ . 

It was argued that the extraordinarily high energy and particle number densities reached in central nuclear 
collisions at RHIC, e 10 - 20GeVfm"^, dN/dy 1000 [||, could lead to rapid (local) thermahzation of 
matter |^ and thus to the creation of the so-called quark-gluon plasma (QGP) a state predicted by 
finite temperature lattice QCD calculations where chiral symmetry is restored and quarks and gluons 
arc dcconfined. 

In order to assess whether this state can actually be formed in an ultrarelativistic nuclear collision one has 
to gain a better understanding of the initial conditions and, at these energies predominantly hard, parton- 
parton scattering processes in the early stage preceding (local) thermodynamical equilibrium. While event 
generators based on individual parton-parton scattering processes Q have been developed, their respective 
predictions for the range of accessible energy and particle number densities differ widely. 

One of the main reasons is the poor understanding of the initial conditions for the nuclear reaction. 
Recently, McLerran and Venugopalan have made considerable progress in a classical approach to construct 
the gluon field at small values of x. Their treatment is somewhat similar to the approach used by Mueller 
for constructing the wave function and gluon structure function of a heavy quarkonium state At small 
X, the nucleonic structure is dominated by gluons, and thus a proper description of gluon dynamics in this 
kinematic region is vital for understanding the initial conditions and the subsequent pre-equilibrium stage 
in nuclear collisions. 
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The McLerran-Venugopalan model considers a very large nucleus moving at ultrarelativistic velocity, 
which consequently appears in the laboratory frame as a "pancake" in the transverse plane. It is assumed 
that due to the large size of the nucleus the (transverse, two-dimensional) color charge density p{x} is large 
(i.e., in a higher-dimensional representation of the color algebra) so that in a certain kinematic region the 
soft gluon field produced by these color charges is effectively classical, and can thus be obtained by solving 
the classical Yang-Mills equations of motion. This field can then be used to compute the gluon distribution 
function. Quantum effects can be implemented as corrections to the classical field. 

The kinematic region for which this approximation is valid is given by the following consideration 
the strong coupling constant as = 5^/47r should be small, therefore the typical gluon transverse momenta 
in the problem should satisfy fcj^ 3> Aqcd- On the other hand, the gluon transverse momenta should be 
sufficiently "soft" , such that the gluons do not resolve individual color charges but couple to the classical 
color charge density. At very high transverse momenta quantum effects become important. Therefore, we 
have to limit ourselves to the region where <^ /z, with /i^ being the average color charge density squared. 
The momentum fraction x of the gluons should be small enough so that the nucleus appears coherent in the 
longitudinal direction. 

To facilitate the inclusion of quantum corrections, the authors of searched for the classical gluon field 
of such a nucleus in the light-cone gauge. The solution of the equations of motion requires to take the 
longitudinal extension of the nucleus into account, i.e., one must not take the nucleus to be infinitely thin 
in the longitudinal direction, as was assumed originally in ||^. The classical field of a single ultrarelativistic 
nucleus is the non-Abelian Weizsacker- Williams field. It was computed in |l^jl^. In the approach pursued in 
PI the nucleus was assumed to be an ensemble of nucleons consisting of point-like (valence) color charges. 
Instead of a smooth two-dimensional color charge density p(x), this quantity is a sum of 6 functions in 
[ pT[ . The limit of applicability of the classical approximation, as well as the structure of the non-Abelian 
Weizsacker- Williams field in terms of Feynman diagrams has been discussed in . That model also allowed 
for an explicit calculation of the average color charge density squared, /x^. 

One of the goals of the approach of McLerran and Venugopalan is to obtain the BFKL equation for 
the structure function of soft gluons, and, if possible, derive corrections to this equation which account for 
nuclear shadowing. Recently, a first step in this direction has been made: the BFKL equation was obtained 
via a renormalization group approach |13| . 

A collision of two ultrarelativistic nuclei of the type advertised in ||^ was considered in Q . The Yang- 
Mills equations were simplified assuming that the solution for the gluon field in the forward light-cone is 
boost-invariant. The equations were then solved perturbatively to first order in the corrections to the Abelian 
solution in the gauge a;_A_|_ + a;+y4_ =0 (where x± — {t±z)/^/2, and A± = (A"±yl^)/V2 are the hght-cone 
components of the gluon field) . 
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In this paper, we also focus on an ultrarelativistic nuclear collision, but we employ the approach of ||ll| to 
describe the nuclei. Each nucleus is moving with the speed of light. The nuclei are taken to be ensembles of 
nucleons, consisting of point-like color charges (valence quarks). The first nucleus is supposed to move in the 
"+" -direction, the second in the "-"-direction, see Fig. 1. In contrast to iQjo), they do have a longitudinal 
extension, i.e., a color charge in nucleus 1 has a fixed a; --component, x^i, which is different for each charge, 
X-i X-k, i k, and similar for nucleus 2, where the charges have fixed x-i- -components, y+j, (of course, 
all charges have different transverse coordinates Xj, as well). 

During the collision, we assume the momenta of the charges to remain unchanged (eikonal trajectories). 
This is certainly justified, since we consider the initial momenta of the charges to be rather large (if not 
infinite). The nuclei just pass through each other and continue their motion along the light-cone (see Fig. |l|). 
As in M, we also solve the classical Yang-Mills equations perturbatively to first order in the corrections to 
the Abelian solution and obtain the classical, radiated gluon field. In contrast to Q, however, we shall work 
exclusively in the covariant (Lorentz) gauge. The advantage of this gauge is that, for the case of a single 
ultrarelativistic nucleus, the classical gluon field is identical to the solution of the corresponding Abelian 
problem Moreover, as will be shown below, for the collision of two nuclei, one can easily relate the 

classical solution to that of a diagrammatic approach in terms of the usual Feynman rules in covariant gauge 
[such rules do not exist in momentum space for the gauge chosen in |^]. 

The outline of the paper is as follows. In Section II we discuss the lowest order solution to the Yang-Mills 
equations. This is just the solution to the corresponding Abelian problem, i.e., the fields generated by the 
color charges in the nuclei are simply superposed and no real gluons are produced. In Section III we find the 
classical field to order g^. This is the first (and lowest) order correction to the Abelian solution. We shall 
derive an explicit expression for the field in coordinate space. We establish the correspondence between this 
classical result and a particular set of Feynman diagrams. This proves that, at this order, the classical field 
is the (major) source of soft gluon production. We shall argue that at higher orders the classical description 
will fail, since already at order non-classical contributions to the gluon field become important. In Section 
IV we calculate energy, number, and multiplicity distributions of the produced gluons. As expected, their 
form is similar to the one found previously inp,|l^. However, the prefactor of our result is different from 
that in while it agrees with the result of [|14[ and [ p9[ . 

Our units are h = c = I, and the metric tensor is g'^'^ = diag(-|-, — , — , — ). Light-cone coordinates are 
defined in the usual way, a± = {a^±a'')/V^, = d/dx±. The notation for transverse vectors is a = (a^^, a^). 



II. CLASSICAL SOLUTION TO LOWEST ORDER IN THE COUPLING CONSTANT 

We consider two nuclei with mass numbers Ai, A2 moving towards each other with ultrarelativistic ve- 
locities, vi^2 — ±1, along the z-axis (cf. Fig. |l|). The nuclei are taken as ensembles of nucleons ||ll[]. In 
order to simplify the color algebra each "nucleon" consists of a quark-antiquark pair. These valence quarks 
and antiquarks are confined inside the nucleons (visualized as spheres of equal radius in the rest frame of 
each nucleus). In order to construct the solution, nucleons inside the nucleus and valence charges inside the 
nucleons are assumed to be "frozen", i.e., they have definite light-cone (and transverse) coordinates which, 
due to our assumption of eikonal trajectories for the individual charges, will not change throughout the 
calculation. We label the coordinates of the quarks in nucleus 1 by x^i, x.^, i — 1,. . . ,Ai, and those of 
nucleus 2 by y+j, y^, j — 1,. . . , A2. Antiquark coordinates follow this notation with an additional prime. 

As in , the nuclei are supposed to be sufficiently "dilute" , such that the distance between the nucleons is 
much larger than the nucleon's size. 

The goal of this section is to solve the classical QCD equations of motion, 

to lowest order in the strong coupling constant (our convention for the covariant derivative is Z?^ = — 
ig[A^, ■]). We shall work exclusively in the covariant gauge {d^Af^ = 0). In this gauge, the equations of 
motion (yj) can be cast into the form 

□A^ = + ig [A^.^d'^A^' + F"^^] , (2) 
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where □ is the d'Alembertian operator. In this form, it is easy to solve the equations perturbatively, as wiU 
be outhned in the foUowing. 

To lowest order in g, i.e., order g, the commutator term on the right-hand side of Eq. (||) does not 
contribute, since the field itself is already of this order [see Eq. (^) below] . To this order and in covariant 
gauge, the classical current can be taken as a sum of the currents for each individual nucleus, as given in 

0: 

Ji'^ =gY,T''{Tn[S{x--x-i)S{x-x,)-Six^~x'_,)6ix-x!,)] , (3a) 

i=l 
A2 

J^y ^gJ^T" (T/) [6{x+ y+,) Six - y.) - (5(x+ - y'^^) 6{x - y^)] , (3b) 

= , (3c) 

where are the generators of SU{Nc), and (Tf) and {Tj") are color matrices which represent the color 
charge of the quark in the color space of nucleon i in the first and nucleon j in the second nucleus, respectively 
[sec [|ll|]. The current (|^) takes into account that antiquarks have the opposite color charge and thus ensures 
color neutrality for each nucleon. 

The classical gluon field satisfying the Yang-Mills equations for a single ultrarelativistic nucleus of our 
type was found in . To lowest order in the coupling constant the solution of Eq. (H) for two nuclei will 

be just a sum of the solutions for single nuclei, since the equations of motion are Abehan {OA'j}^ = jji'^), 
and thus linear. Therefore, as one readily confirms by an explicit calculation, the solution of the Yang-Mills 
equations to order g is [|ll| : 



Ai 



4' --|^E^"(^") - ^^r) Hk~ - Six- - x'_,) \ni\x- x'^\X)] , (4a) 



i=l 
A2 



A^^^ = -fJ2^^f^) [Six+-y+,) H\x-y^\X)-S{x+-y'^^) ln(|ai - ^;|A)J , (4b) 



2tt 

^(1) = , (4c) 

where A enters as an infrared cut-off. In a sense, it acts as a gauge parameter that sets the scale of the gauge 
potential. The associated field strength tensor is independent of A, 

(5a) 

a:_,)-^^-S{x--x'_,)j^^4^] , (5b) 



^1) 


-0 , 






g 

2tt 


Ai 

E 




_ 9 
27r 


A2 

E 




= , 





\x — x-^^ 



x — y. 2L^ y' ■ ' 



(5c) 



(5d) 



We note that the field strength is zero in the forward light-cone. This is of course reasonable, because to 
order g there are no interactions between the point-like charges constituting the nuclei. Therefore, no real 
gluons are produced. 

III. GLUON FIELD TO NEXT-TO-LOWEST ORDER 

A. Formal solution of the equations of motion 

In this section we compute the solution of the Yang-Mills equations to order g^. As we shall see, the 
equations of motion (j^) are linear to each order, and we therefore focus first, for the sake of simplicity, on 
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the case of a collision of two single point-like color charges, for instance the quark from nucleon i in nucleus 1 
and the quark from nucleon j in nucleus 2. The generalization to the nuclear collision is then straightforward: 
the solution is a simple superposition of the solutions emerging from each individual collision (i.e., a sum over 
i and j and over the respective quark-quark, quark-antiquark, antiquark-quark, and antiquark-antiquark 
scatterings). 

To order g'^ the equations of motion (|^) read: 

where ^{^^^ and J^i^' are the contributions to the giuon field and the fermionic current to order g'^. In order to 

solve these equations we have to first determine Jfj. . The most simple approach is to exploit (non-Abelian) 
current conservation J^^ = 0. One obtains: 



3 

J(3)^' = tg [A« , + ^g [aL') , Jj'^] = I- r"' T'^ (T^) (T;) 

ZTT 



X S{x^-x^,)6{x+-y+j) S{x-y)-6{x-x,) ln(|xj - y |A) . (7) 



The charges are assumed to be recoilless and follow eikonal trajectories. Therefore, their momenta do not 
change in the interaction and the transverse component of the fermionic current is zero. The "+" and 
"-"-components will still be d functions on the light cone and in transverse direction, as was the case at 
order g. The only effect of the collision on the valence charges is a "rotation" of their color, as soon as a 
charge "hits" the field of the other charge at the collision point. This consideration leads us to the conclusion 
jf'' ~ (5(x_ —x^i) 9{x+ — yj^j) 5{x — Xj) and J^^ ^ 0{x- —x^i) 5{x+ — y+j) Six — y.). The correct coefficients 
are found from Eq. (Q): 

Jf = -^^r"' T'^ (T^) (T^) H^- - ^-i) - V+j) HSL - Hi) ln(k, ~ y^. |A) , (8a) 
jiZ) ^ |!/«&CT''(Tf)(T;^)0(a:_ -a;_.)<5(x+- ?;+,)% -yp ln(|a:, - ^^.|A) , (8b) 
J(3) = . (8c) 

It can be shown that this current is consistent with the eikonal scattering limit of either the lowest order 
QCD diagrams for gluon radiation from two colliding color charges (cf. also Sect. IIIB below) or Wong's 
equations jl^ for the collision of two classical color charges. 

For one charge, say i, "hitting" an ensemble of charges j, the color of charge i "rotates" each time it hits 
the field of one of the charges j in the ensemble. The resulting current is simply the sum over j of the terms 
on the right-hand side of (ph . Similarly, the collision of an ensemble with an ensemble simply adds another 
summation over i. The actual nuclear collision in our approach is only slightly more complex in that one 
has to account for the presence of antiquarks as well. 

The equations of motion for the next-to-lowest order gluon field are now obtained by inserting the lowest- 
order results together with the current (||) into the right-hand side of Eq. (^): 

□ Af'' - ^ (T^) {f^) (- 2^ ln(|x, - yj\) 5{x^ - x^.,) 9{x+ - y+,) 5{x - x^ (9a) 

+ ^'^^{\x-Xi\X)\n{\x-y^\\)d+5{x--x^i)6{x+-y+j)^ , 

□ = (1^ (^') ( - ^,1^) ^(^- - ^-^) '^(^+ - y+i) - V^) (9b) 

- ln(|x - Xj|A) \n{\x-y ]^X)5{x^ - x-t)d-5{x+ - y+j] 

□ A^'^"" = ^ r"' (Th {T^) Six- - x_.) Six+ ~ y+,) (9c) 

X \ H\x-y\X) - -Inds-xJA) ~ 

\ J \X XA 




5 



Equations are linear differential equations, which justifies why we were able to focus on one single 
collision between two valence quarks first and later obtain the complete solution by summing over all possible 
collisions between color charges. (It was explained above how the corresponding sum over i and j and over 
quarks and antiquarks appears in the fermionic current; for the commutator term, its presence is obvious.) 
The linearity of (^) also allows us to compute the solution simply by the method of Green functions. Since 
the classical solution obeys causality, we have to use the retarded Green function: 

(a;) = J d^x'Grix - x') J^(x'), (10) 

where = + tg [A^^^, 9,4'^ + Pij}] is given explicitly by the right-hand side of Eq. (^. The retarded 
Green function reads in coordinate and momentum space : 

Gr{x)^^0{t)S{x^) , Grik) = - -g^^-j- . (11) 

Formulae (^ ~ ^I]) provide the classical gluon field to order g'^. Note that the above perturbative solution 
scheme renders the equations of motion linear at each successive order in g. In principle, one can therefore 
use the method of Green functions to construct the classical solution to arbitrary order in g. We shall argue 
below, however, that already at order g^ quantum effects become important and the classical approach breaks 
down. Before we compute the solution to order g^ explicitly, let us draw a connection to the perturbative 
solution via Feynman diagrams. 



B. Connection to Feynman diagrams 



Let us write the right-hand side of Eq. (hOh in momentum representation 



4)«(. 



(27r)4 fc2 + ieko 



Juik) 



where 

Jl{k) 
l\k) 



9 fabc (rpb\ /tSc\ + ■) 



(fqe 



-iq-(x.-y .) 



{k-qf 
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(2^)2 

„3 



ik_ 



ik-qy 



i2ny 



a qe - 



i{2q-'k) 
q'{k-qy 



In order to derive this, we have made repeated use of 



d^g ■ 1 1 



(27r)2 



27r 



(12) 

(13a) 
(13b) 
(13c) 

(14) 



and, for (13c), of the transverse gradient of this equation. 

The diagrams giving the classical gluon field at order g'^ are shown in Fig. ||. The cross at the end of 
the gluon line denotes the space-time point x where we measure the field. The upper quark line corresponds 
to the first charge (fixed coordinates x_j), the lower one corresponds to the second charge (at fixed 
Vj)- Therefore, the momentum on the upper line has a large -component, the momentum on the 

lower line has a large "-"-component. The gluon field to order g'^ also includes graphs where the two quarks 
do not interact. These diagrams arc not shown in Fig. ^, since they are not part of the classical gluon field 
and do not contribute to gluon production (they vanish once we take the emitted gluon line to be on-shell). 
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We take the gluon-fermion vertices in the eikonal approximation. In the standard calculation of the 
diagrams, the emitted gluon line corresponds to a gluon (Feynman) propagator — i/(fc^ + ie) times a phase 
^-ik-x ^gjj^(2g TffQ shall ultimately transform the diagrams into coordinate space). This term is common 
to all diagrams in Fig. ^ However, due to the regularization of the Feynman propagator, all diagrams 
corresponding to gluon absorption (instead of emission) are automatically included in this calculation. That 
means, the usual Feynman diagrams yield an acausal result. In order to establish correspondence to the 
classical result, where there is only gluon emission, we replace the Feynman propagator by the retarded 
propagator —i/ {k^ + ieko) for the gluon field measured at x. 

Similarly, due to the Feynman regularization of the fermion propagator, all diagrams B — E contain 
contributions where the gluon is emitted prior to the one-gluon exchange. To ensure causality, we have 
to use a retarded fermion propagator when calculating the diagrams B and D and an advanced fermion 
propagator in diagrams C and E which renders the emission of the gluon causal. 

We have the freedom to change the regularization of propagators. Different choices of regularization do 
not influence the physics. The difference between the retarded (or advanced) propagators and the usual 
Feynman propagator is just a d function of the square of the four-momentum of the internal line. Therefore, 
for graphs B — E this difference is proportional to those parts of the diagrams where the fermion line is 
on-shell. But these parts do not contribute to real gluon production, since once we put the outgoing gluon 
on-shell, they vanish. The regularization of the internal gluon lines of the diagrams in Fig. |2| turns out to 
be of no importance for the actual calculation. 

After clarifying the regularization of the propagators, we compute the diagrams according to the usual 
Feynman rules in covariant gauge. Let us denote the pth component of the gluon field from a diagram X 
by X\p, where X = A, B,C, D, or E. Diagram A is non-zero for all values of p, while for B and C only 
the "-f "-component and for D and E only the "-"-component are non- vanishing. After a lengthy, but 
straightforward calculation we compare with Eqs. (|2|) and ( [l3| ) to obtain the identities 

Af'^(a;) = (A + S + C)|,^+ , (15a) 

A'^^^'^ix) = iA + D + E)\p=- , (15b) 

A(3)-(x) = A\p^^ . (15c) 

We see that the calculation of the diagrams yields exactly the gluon field ( p^ obtained from the classical 
solution of the Yang-Mills equations. That proves the correspondence of the classical field to the diagrams 
in Fig. ^. The diagram A arises from the commutator term on the right-hand side of Eq. (^), whereas graphs 
B — E arise from the fermionic current. 

The order is the limit of applicability of the classical approach to the problem of gluon production. 
At order one can construct diagrams which contribute to gluon production, but cannot be obtained 
classically. An example of this type is shown in Fig. 0. 




FIG. 3. Graphs contributing to the gluon production at order g^. 
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The two-gluon exchange contribution (Fig. ^ calculated using the traditional Feynman regularization of 
the propagators cannot be obtained from the classical equations of motion. The classical approach would 
give the part of the diagram where the two exchanged gluons are in a color singlet state. That corresponds 
to the internal fermion lines between the exchanged gluons being on-shell. This is what one would obtain 
by iterating the procedure for the determination of the classical field outlined above to higher orders in the 
coupling. The color octet combination of the two gluons, being the lowest order correction to the gluon's 
Regge trajectory, is a pure "quantum" effect. A different regularization of the propagators cannot make this 
contribution "classical" . One can see that by a direct calculation of the two-gluon exchange diagrams in 
the octet channel. The resulting contribution is proportional to In s, where s is the center of mass energy 
of the colliding quarks [see and references mentioned there]. That is, it depends on the longitudinal 
momenta of the quarks. However, the classical calculation is just an iteration of the one-gluon exchange 
diagram and thus not able to provide such a logarithmic longitudinal dependence. Therefore, the diagrams 
in Fig. H are not classical. (After all, it would be quite surprising if the gluon's Regge trajectory turns out 
to be a classical concept.) 



C. Classical gluon field in coordinate space 

Let us now explicitly compute the gluon field ( p^ ) in coordinate space. We decompose the "+" -component 
of the field as 



A\' = fli + 02 , 



(16) 



where ai is straightforwardly computed: 



ai 



9 fabc trr,b\ 



9fabc (y.) (yc) L 



2(27r)- ■' ■ ~ X _ 
The second term can be written in the form 



X- - .T_j) 9[x+- y+j 



2{x- - x-i) 



H\x,-yJ\) . (17) 



"2 = ^ {T^) if^) J d^x' Gr{x - x') H\x! - x^X) ln(|x' - y^.|A) d'+5{x'_ - S{x'+ - y+,) 



9 jabc (yc) Q^^g^^_ _ ^(^^ _ y^^^j^ 



2(27r) 



with 



_ j' Cpq(Pl -iq-{x-y,)-il-(x-x.) 1 



(18) 



(19) 



where r — ^J2(x^ — x_i)(x+ — y+j). The explicit evaluation of the integral J is referred to Appendix A. 
The final result is: 



^-ln(C>A) ln(7;>A) + i 



Li2 e 



£,>V> 



+ Lis e 



(20) 



where C>(<) = max(min)(|a; — xj, r) , ??>(<) = max(min)(|a; ~ Vjl^ a is the angle between x — and 

x — y.. A plays the same role as in Eq. (^). Li2(z) is the dilogarithm (also known as Spence's function). 
With the final expression for 02 we obtain: 



9{x- - x-i)9 lx+ - y+j - -j^—^ 



ln(|£, - yJA) 



5{x^ ~ x_i)e{x+ -y+j) ln(|x-£jA) ln(|x-y.|A) -e{x_ -x_i)e{x+ - y+j) 9+[ln(C> A) ln(77>A)] 



+ - e{x^ - x^i)e{xj^ - y+j) {d+ Inr) ln(l - 2rcosQ; + r^) 



(21) 
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where r = {£.<V<)/{^>V>)- 

The integrations for A^^" are done similarly with the result: 



1 



\x - I 



2(a;+ - y+j) 



ln(k, - y .|A) 



6{x+ ~y+j)0{x_ -x-i) ln(|^-x,|A) \n{\x-y\X) - 9{x_ -x_i)9{x+ - y+j) (9_ [ln(^> A) ln(77>A)] 



6{x+ — y+j) (9- Inr) ln(l — 2r cosa + r^) 



(22) 



For A^^)" we find 



^^'^'^(^) = {Tt) (f/) 9{x^ - X-,) 9{x+ y+j) [(V. - V,) J] 



2(27r) 



--[(V, - V,) Inr] ln(l - 2r cos a + r^) - - [(V, - V,>] In ( ^— -— 



1 - re*" 



(23) 



As expected, the solution to order is causal, i.e., it is non- vanishing only in the forward light-cone. 

To obtain the complete solution for the nucleus-nucleus collision we simply sum over all charges in the 
nuclei, i.e., as discussed above, we sum over all nucleons i and j and over quarks and antiquaries inside the 
nucleons. If we denote the solution of the scattering problem of two single charges with coordinates Xi and 
yj found above by A''^^°'{x,Xi,yj) we obtain 



Ai,A-2 



(24) 



The antiquark coordinates are marked with a prime. The relative signs emerge from the fact that antiquarks 
have the same color charge as quarks, but with opposite sign. 



IV. THE RADIATED FIELD ENERGY, NUMBER, AND MULTIPLICITY DISTRIBUTIONS 



In order to determine the radiated field energy we start from Eq. (y_2|), and, for simplicity, discuss the case 
of a single collision first. Note that a part of the solution corresponds just to a change of the field carried 
by the charge due to the collision (the color "rotation"), and not to the radiated gluon field. That part is 
most easily isolated by a contour integration in the complex /cQ-plane in Eq. (12), it arises from the pole 
k± = —it in Eq. (p^): 



l(3)a 

+ charge 



{x) = 



fate (yc) 1^(1^^ _ ^1^) 1^(1^ _ ^^1^) _ ^(^^ _ y^^^ ^ 



(27r)2 ^ J' — 



A 



(3)a 
' charge 



(x) 



(27r)2 



r'^iTl) (rn ln(|x,-y lA) H\x-y\X)0{x^-x^,)Six+-y+,), 



(-)=o 



(25a) 

(25b) 
(25c) 



The radiated gluon field arises from the poles of the retarded propagator in Eq. (|l^), i.e., it corresponds to 
on-shell giuons, as one would expect: 



A 



(3)a 
fi rad V 



(a;) = 9{t-t,j) / dk i J°(w,k) e 



(26) 



where tij = (x_i + y+j)/V2 is the time when the collision happens and dk = d'^'k/[{2n)^2Lij], with 
p6| . Obviously, this field vanishes prior to the collision. 
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The (stationary part of the) radiated field energy is ||T^ 



n= dkuj 



J iuj,k)-J (w,k) - J^(tj,k)>*(cj,k) - J^(cj,k) J^*(w,k) 



(27) 



In the case of a nuclear collision J^{k), Eq. (p^, becomes more complicated: there is an additional sum over 
nucleons i and j and over quarks and antiquarks inside these nucleons. Inserting the resulting expression, 
for times after the last parton-parton collision we arrive at 



abc rade 



Ai A2 



dk Lo 



1 



(28) 



X E Y.(^'^(^D(T^nfnr{k+,q^;x,)V{k^,k~q^-,y,)V^k+,q^;Xk)V*{k_,k-q^,yi) 
i,k=lj,l=l 



where now k± = (a; ± fc^ )/ v 2 and 



(29) 



In order to achieve color neutrality in the initial state, we have to average over all possible color orientations 
in the color space of the individual nucleons. With tr [ (Tj^) (f/) (T^) {ff) ] = 6^k S'"^Sji 5^=74, and f'"' = 
Nc{N^ - 1) this yields 



4 tr[H] = 



(2^)4 Nc 



where 



Ai A-i 
Vik,q^,q^) = Y,'P{K,g{,x,)V*ik+,q^;x,)J2'Pik-,k-q^-,y,)r*{k_,k-q^-,y,) . (31) 
i=l j=l 

For further evaluation we introduce the center-of-mass coordinates of nucleon i, (X_i,X,j), and nucleon j, 
(Y^j,Y_j), and relative coordinates Aa;_i, Axj, Ay_(-j, Ay ., such that 



X-i = 



Aa;- 



X. 



Ay 



+1 



2 ' — ' 2 ' 2 ' -J -J 2 



(32a) 
(32b) 



Here we assumed that the positions of quark and antiquark in a nucleon are symmetric with respect to the 
center of the nucleon. This assumption is not crucial. The calculations can also be done for the general case 
of arbitrary positions of quarks and antiquarks in the nucleons. For the "symmetric" case the phase factor 
becomes 



V{Kq^,q^)=Y,e-^^l^-^.>^^ 



i=l 
A2 

X N e -1 -2' ■ 



-i(q -q )-A£C./2 ik+Ax-i-i(n +n )-Ax./2 , 



p — 1 — 2 —J — g —J —1 —2 —J -|_ , 



(33) 



We shall now first average V over the longitudinal positions of quarks and antiquarks inside the nucleons. 
The nuclei are highly Lorentz-contracted in the longitudinal direction. To simplify the subsequent discussion, 
we take them and the nucleons inside to be cylindrical. Each nucleon is then a cylinder of radius a and 
length 2a/ J, oriented along the z-axis, where 7 is the Lorentz-factor in the CM frame of the collision. (The 
longitudinal extension of the nucleons is not important, it will drop out in the following anyway.) With 
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\/2a J-a/^-i \ 2 y k+\/2a/^ 

and an analogous relation for the average over Ay_|-j , the longitudinal momentum dependence vanishes from 
the phase factor. 

Let us now average over the transverse positions X_^ , Y_j of the nucleons inside a nuclear transverse area. For 
a cylindrical nucleus, the transverse area which we average over is independent of the longitudinal position 
of the individual nucleon and equal to irR^^^y i?i(2) being the transverse radius of nucleus 1 (2). Let us take 
nucleus 1 to be the larger of the two nuclei. Since Ri is by far larger than the (inverse) momentum scales 
we are interested in, we may take to good approximation i?i — )■ oo, and obtain 

,35, 

With this result, the F^-average (over the transverse positions of nucleons in nucleus 2) is trivial. The 
assumption of infinite transverse nuclear size greatly simplifies the subsequent discussion. However, we are 
then no longer able to study collisions at finite impact parameter. 

Finally, we average over the transverse dimension of the nucleons. That integral is formally the same as 
(p5|), only that Ri is replaced by a and the momentum dependences of the various terms in V are different. 
Now, however, 1/a is not small on the momentum scale of interest. In fact, the ^"iiitegrations in (|3^) are 

logarithmically divergent, such that the treatment of the small momentum region is of some importance. 
As we shall see, treating the average over the transverse dimensions of the nucleons correctly leads to an 
infrared regularization of these divergences via nucleonic form factors. The physical interpretation is that 
on large (spatial) scales individual nucleons appear colorless and do not emit gluons. 

The average over the transverse dimension of a nucleon involves typically integrals of the type ||l8| 

na^ J \ 2 J \q\a 
For the average phase factor we thus obtain 

47r / 2Ji(2|gla)\ / 2 Ji{2\k - q\a)\ 

Our final result for the invariant distribution of the average energy is 

d{n) g6 - 1 4A,A2 u f F{\q\a) F(|fc-g|a) 

dyd^k ^(2^)6 TTRj J ^ f (k-q)^ ' ^ ' 

with the (longitudinal) gluon rapidity y = ln[A:+/fc_]/2 and the "form factor" F{x) = 1 — 2 Ji(2x)/2x. For 
small q or k — q, F regularizes the infrared divergence of the integral. If we do not assume that the quark 
and antiquark are symmetric with respect to the nucleon's center and allow them to be anywhere in the 
nucleon, the "form factor" becomes F{x) ~ 1 — [2 Ji{x) / x]"^ . 

We now compare our result to that of [||. In that work, the number distribution of radiated gluons was 
related to the average energy distribution via 



dJV _ 1 d{H) 



dycPk to dy d^k 



(39) 



As in 1^, we consider a central collision of equal nuclei, Ai ^ A2 ~ A. According to ijlTj, the average 
transverse color charge density is given by /i^ = CpA/ {NcTtR'^) (for our case of a "cylindrical" nucleus), 
where Cp = (N^ — l)/{2 Nc). Then, with the transverse area S±_ = Tri?^, 

dM ^ TV 25V ,^t2 .^ 1 /■ F{\k~g\a) 



dyd^k Cl (27r)37r " ' ^ J (27r)2 q^ {k~q) 



S^VW7it^Nc{Nt~l) -2 7:^ r J ■ (40) 
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This has (apart from the form factors) the same momentum dependence as the result of § [Eq. (49)]. The 
prefactor, though different from agrees with recent results obtained by Gyulassy and McLerran | |l9|] . 
Note that due to the form factors, our result behaves like 1/fc^ for smah |fc| ~ 1/a. The residual infrared 
divergence stems from our approximation of an infinitely large nucleus in Eq. (Bq) . This divergence is actually 
cut off by the finite (inverse) size of the nucleus (i.e., by the nuclear form factor). In this range of momenta, 
however, our treatment ceases to be valid anyway, since then |fc| ~ 1/R <ti 1/a ^ Aqcd, while one of our 
initial assumptions was |fc| ^ Aqcd. 

The gluon number distribution ( pO[ ) is not directly accessible in the experiment. Experimentally, one 
measures the (invariant) differential cross section for detecting a radiated gluon, da-^i^^l dy dPk, divided by 
the total cross section of the collision. Equivalently, one can divide the gluon number distribution ( pO| ) 
by the total number of scattering events with one-gluon exchange between the color charges in a (central) 
A + A-coUision. This number can be obtained quite analogously to the above derivation of the number 
distribution for radiated gluons. Instead of Eq. (2^) we start with the square of the one-gluon exchange 
amplitude. This amplitude for two charges (whose coordinates are specified as in Sect. Ill) is 



-f= e - -J 



(41) 



Similar to the above, after averaging the modulus squared of this amplitude over the positions of the quarks 
and antiquarks inside the nucleons and nucleons in the nucleus, as well as over quark colors in the initial 
state, we obtain 



2 Cf 4^2 

4 as 



d^q 



2Nc 7ri?2 

The final answer for the gluon multiplicity distribution is 



Fi\q\a) 



(42) 



dn 



1 



dAf 



dyd^k Mtot dyd^k k 



Ncas 

2 1,2 



d^q 



Fi\q\a) F{\k^q\a) 



q- 



{k-qf 



d^q 



(43) 



Expanding the form factor (in the "symmetric" case) for small a, F{x) ~ + a;^/3) for a; — > 0, and 

choosing = 24/0^, this result coincides with that of [|4j [Eq. (21)], since Ca = Nc = 3. 

For the case of a nucleon in which the positions of quark and antiquark are symmetric with respect to its 
center one can explicitly perform the integration in Eq. (H3). The result reads: 



dn 
dy d'^k 



N^as 12 f 1 + (kaf 



'k' 



{kflY 



(1- Jo(2|fc|a)) 



(fca) 



[kaY 



The details of the integration are given in Appendix B. For k 
becomes 5/12, as expected from (Esf). 



J2{2\k\a) + 2F3(l,l;2,2,2;-fcV)j . (44) 
0, the expression in large parentheses 



V. CONCLUSIONS 



In this paper we have solved the classical Yang-Mills equations for a collision of ultrarelativistic nuclei. 
Our solution assumes the nuclei to have the form used in ||ll[ and the color charges to follow eikonal 
trajectories during the collision. The solution is constructed perturbatively in covariant (Lorentz) gauge to 
lowest order (Abelian limit, order g) and next-to-lowest order {g^). Via Feynman diagrams we have clarified 
the connection between the classical and quantum solution to order g^ . 

We have established a limit for the classical description of gluon production. From the discussion of the 
diagrams in Sect. IIIB follows that the limit for the field is one gluon per nucleon, unless the second gluon 
is the outgoing gluon. That means that in the diagrammatic formulation of the problem there should be no 
more than one gluon leaving each nucleon. In nuclear collisions, different from the case of a single nucleus 
[^|, we cannot require the final states of the nucleons to be color singlets. Therefore, here we cannot apply 
color averaging which cancels quantum corrections [as described in [p2[]] and puts internal fermion lines of 



12 



graphs like in Fig. y on shell, allowing for a classical description. That way the limit of the classical approach 
is reduced to one gluon per nucleon. 

We have given explicit expressions for the gluon field in coordinate space [Eqs. (|2| - [2p)] and calculated 
the energy, number, and multiplicity distributions of radiated gluons to order [Eqs. (|38|, ^]. Our 
result for the number distribution agrees with |]l9| and the form factors of the individual nucleons can be 
cast into a form such that the multiplicity distribution agrees with the result of p^ . 

The resulting gluon number and multiplicity distributions are boost-invariant, a property which was a 
priori assumed in Boost-invariance is also assumed in popular hydrodynamic models to describe 
the evolution of ultrarelativistic nuclear collisions after (local) thermalization is established. To answer the 
question whether thermalization actually happens and what the respective time scales are, it is planned to 
utilize the solution found here to study screening and damping in the radiation-produced gluonic medium 
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APPENDIX A: 

In this appendix we explicitly calculate the integral 

^ = j^ P^i + ^1^)- (^^) 

Using the addition theorem for Bessel functions [(j) is the angle between q and /) 

oc 

M\i+l^T)^ ^ (-l)V'=^Jfe(MT)Jfc(|l|r) (A2) 

fc— — oo 

we obtain: 

^= E (-1)' ih^ ir-, e^'^M\q\T)M\l\r). (A3) 



k— — c 



If 01 is the angle between I and x — x^^, and (f)2 is the angle between q and x — y.^ we can express <f> as 

= 02 - 01 + a, (A4) 

with a being the angle between x — y ^ and x — x^. All angles are taken clockwise. Integrating over 0i and 
02 we obtain: 

jMq\3L-y^\)Jkiqr)^ (^J^^jJk{l\x-x,\)MlT)^ . (A5) 
For non-zero k we can use the formula [see [ p2[ ] 

'^^■Jk{l\x-x^\)MlT)^^(^) , k>0, (A6) 
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where C>(<) = ni3'X(min)(|a; — xj, r), to obtain the series: 



k=l 



with '7>(<) = max(min)(|x ~ y ],t). Using the definition of the dilogarithm, 



Li2(^) = X]f2 



fc=i 



we can rewrite this series as 



Lio e 



Lio e 



For the (fc — 0)-term in Eq. (|A5|) one can show that since 



f 



dl 



Jo{1\x~x^)Jo{It) = - hi 



with 7 being the Euler constant, this term can be rewritten as 

ln(^>A)ln(r;>A) , 



(A7) 



(A8) 



(A9) 



(AlO) 



(All) 



where the cutoff A includes the previous cutoff /i and all other numerical prefactors. We finally obtain the 
answer: 



J = ln(e>A)ln(r7>A) + i 



Li2 e 



A<V< 
£,>V> 



+ Lis e 



£.>v> 



(A12) 



APPENDIX B: 



The goal of this appendix is to perform the following integration: 



I = 



1 - 



.h{2\q\a) 



1 - 



Ji{2\k-q\a) 



g2(fc-g-)2 \q\a J \' |fc-g|a 

Defining g = j^j and p= |fc — g] we rewrite the integral as 

r J{p,q)dpdq / _ Ji{2qa) \ / _ M2pa) \ 
J q^P^ V qa ) \ pa ) ' 

where the Jacobian J(p, q) is given by [see [p|,p2|] 



J{p,q) = 2npq / bdb Jo{bk) Jo{bq) Jo{bp) , 
Jo 



with k = \k\. Inserting Eq. ( B3 ) into Eq. ( JB2| ) we obtain: 

I=2n HbdbMbk) ir^Mbq) (l-'Il^) 

Jo Uo 1 V go y 



(Bl) 



(B2) 



(B3) 



(B4) 



The integration over p in Eq. (B2) is identical to the integration over q, which allowed us to square the 
g-integral in Eq. (B4). Performing the integral in the square brackets we get [see p|,p2[] 
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jMbq) [ 1 



Ji{2qa) 



qa 



which yields for the original integral 

r2a 



1= 27r / bdbJoibk) 
Jo 



e{2a-b) 



1 

2 ^ 8a2 



(B5) 



27r(2a)2 / t dt jQ{2akt) [Xnt + -{I - t^) 



(B6) 



with t = 6/2a. The integral in Eq. (B6) can be calculated by expanding the Bessel function in a power series, 
performing the integration in each term, and finally resumming the series. The result is 



X = 27ra'^ 



^ ^ ^^"^^^ (1 - Jo{2ka)) + M2ka) + 2^3(1, 1; 2, 2, 2; -fcV) 



(fca)4 



(fca)4 



This expression is used to obtain Eq. ([44|). Taking the k ^ Q limit of formula 
appearing in the number of elastic scattering events [Eq. 



d^q 



(g2 



=T 1 



Ji(2|g|a) 



— 7ra 
6 



(B7) 

we obtain the integral 
(B8) 
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